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Abstract 


Within the framework of perturbative Quantum Chromodynamics we derive transverse momen¬ 
tum dependent distributions describing both current and target fragmentation in semi-inclusive 
Deep Inelastic Scattering. We present, to leading logarithmic accuracy, the corresponding cross- 
sections describing final state hadrons on the whole phase space. Phenomenological implications 
and further developments are briefly discussed. 


I. Introduction 


The predictivity of perturbative Quantum Chromodynamics (pQCD) relies upon the factor¬ 
ization of hadronic cross-sections into perturbative process dependent coefficient functions, uni¬ 
versal perturbative evolution equations and non-pertnrbative process independent densities [1]. 
In presence of a hard scale, which justifies the perturbative approach due to asympotic free¬ 
dom, the scale dependence of densities is predicted through renormalization group equations, 
i.e. Altarelli-Parisi (AP) evolution equations [2]. Inclusive deep inelastic lepton-hadron scattering 
(DIS) has been investigated since a long time and in great detail. The evaluations of the splitting 
functions has been recently performed np to O(a^) in perturbation theory [3], further constrain¬ 
ing non-pertnrbative dynamics. The evolution of initial state partons in terms of longitudinal 
and transverse momenta has been also considered within pQCD in Refs. [4-6]. In semi-inclusive 
Deep Inelastic Scattering (SIDIS) processes, at variance with inclusive DIS, one hadron is de¬ 
tected in the final state, l + P —> l + h + X. In this case, on the contrary, an equally accurate 
theoretical description has not yet been developed. The additional hadronic degrees of freedom 
require a more detailed description of parton dynamics. Within the usual pQCD-improved par- 
ton model approach to SIDIS [7,8], one deals only with current fragmented hadrons. However 
another distinct issue enters the perturbative description. Both the struck parton and spectators 
do evolve according to the hard scale governing the process [9]. As a result also target fragmen¬ 
tation has to be included to describe final state hadrons. Evolution is predictable in terms of 
new functions dubbed fracture functions, whose factorization can been proven in Ref. [10]. An 
explicit evalution of the single-particle cross section at one loop has been given in Ref. [11], 

In the pQCD-improved parton model detected hadrons are espected to have a sizeable trans¬ 
verse momentum Ph±, as a result of perturbative evolution in terms of hard partons emission. 
We will show that it is possible to reformulate evolution equations in order to include transverse 
degrees of freedom since such a dependence is fully predictable within pQCD. 

The aim of this work is twofold. First we derive transverse momentum dependent (TMD) 
evolution equations which enter SIDIS cross-sections in the current fragmentation region. We 
then extend this treatment to distributions in the target fragmentation region via fracture func¬ 
tions. As a result of these generalization, the combined SIDIS cross section is presented descibing 
hadron production on the whole phase space. 


II. Transverse evolution equations and kinematics 


Evolution equations to leading logarithmic accuracy (LLA) resum contributions due to collinear 
partons emission. In the time-like case the evolution equation is: 
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where the fragmentation function V^(zh) represents the probability to find, at a scale Q 2 , a given 
hadron h with momentum fraction Zh of its parent parton i. P{u ) are the time-like splitting func¬ 
tions which give the probabilities to find a parton of type j inside a parton of type % and can be 
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expressed as a power series in the strong running coupling a s (Q 2 ) . Ordinary evolution equations, 
eq. 0 , contain only longitudinal degrees of freedom of partons inside hadrons although, at each 
branching, the emitting parton acquires a transverse momentum relative to its initial direction. 
Transverse momentum dependent (TMD) evolution equations were first derived in Ref. [12] for 
fragmentation functions in the time-like region: 

Q* aV ‘ (Z ^’ P±) = 

■PlZ 5(u(l - u)Q 2 - q\)T*(—,Q 2 ,p± - —q±). (2) 

Single particle distributions V^(zh) in eq.© give the probability to find, at a given scale Q 2 , 
the hadron h in the parent parton i with longitudinal momentum fraction Zh and transverse 
momentum p± with respect to it. The P(u) splitting functions are the ordinary splitting functions 
as in eq.© and flavour indices are understood as in the inclusive case. Let us discuss the 



Figure 1: a parton with momentum k emerges from a hard process with virtual mass k 2 = Q 2 > 0 and 
then evolves into a quasi-real parton with transverse momentum p±. The vertex is associated with the 
time-like splitting functions P(u). The small blob symbolizes resummation of ladder diagrams in LLA. 


kinematics structure of eq. ©• Consider an outgoing parton k emerging from an hard collision 
with virtuality k 2 = Q 2 > 0, assigned to have zero transverse momentum (fc_|_ = 0j_) and 
unitary longitudinal momentum fraction, as displayed in Fig.(l). ft subsequently branches, with 
probability P(u), into a couple of partons, q and q', with a fractional momentum u and 1 — u 
of k respectively and transverse momentum q± = —q'^ relative to k. The parton q then non- 
perturbatively hadronizes generating the final hadron h with a fractional momentum Zh and 
transverse momentum p± and p± relative to k and q respectively. We thus derive the following 
constraints: 


Rl = P±-—q±, (3) 

u 

ql = u(l-u)Q 2 . (4) 

Eq. @ takes into account the Lorentz boost of transverse momenta p± from the g-frame to the 
frame. The second equation follows by imposing the virtuality flow of time-like branching. 
These relations directly enter eq.©, respectively as argument of the distribution T>j and of the 


3 










Figure 2: a quasi-real parton with momentum p perturbatively evolves towards the hard vertex acquir¬ 
ing a virtual mass p 2 = Q 2 < 0. The vertex is associated with the space-like splitting functions P(u). 
The small blob symbolizes resummation of ladder diagrams in LLA. 


invariant-mass conserving 5-function. The unintegrated densities fulfil the normalization: 

f d 2 Pl _Vl(z hl Q 2 , P± ) = V?(z h ,Q 2 ), (5) 


since the boost in eq. © is linear in the transverse variables, i.e. the Jacobian is: 
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This property garantees that we can recover inclusive distributions, eq. ©, starting from less 
inclusive ones. The opposite statement is not valid since eq. © contains new physical information. 

In order to obtain a complete description of semi-inclusive cross-sections we need the space¬ 
like version of eq. ©• On a general ground we may expect that the infrared structure of space-like 
evolution equations is the same as the time-like one since it depends only upon the dynamics of 
the underlying gauge theory, the only changes being in the kinematics. In analogy to the time¬ 
like case we consider now a initial state parton p in a incoming proton P which undergoes a hard 
collision, the reference frame being aligned along the incoming proton axis, as in Fig. (2). The 
boost of transverse momentum and the invariant mass-conserving constraint are in this case: 

fcx - «X (7) 
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We thus generalize eq.© to the space-like case: 
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The constructive proof of eq.© will be given elsewhere [27]. The unintegrated distributions fulfil 
a condition analogous to the one in eq. ©, i.e. : 


d 2 k ± F P (x B , Q 2 , k ± ) = Pp(x B , Q 2 ) 


( 10 ) 
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We note that the inclusion of transverse momenta do not affect longitudinal degrees of freedom 
since partons always degrade their fractional momenta in the perturbative branching processes. 
The different Lorentz structure in the transverse arguments of the parton distribution functions 
T' l p arises from the different structure of the Bethe-Salpeter ladder used to derive the evolution 
equations [27]. 


III. Transverse momenta in target fragmentation region 


Current and target hadron production mechanisms cannot be separated since hadrons pro¬ 
duced by current fragmentation may go in to the target remnant direction and vice versa. In 
these configurations new infrared singularities appear which cannot be reabsorbed through the 
standard renormalization procedure into parton distribution functions and fragmentation func¬ 
tions. It has been shown [9-11] that the cross-section can be renormalized by introducing new 
non-perturbative objects, fracture functions indicated by A 4 Ph (x, z, Q 2 ). These functions express 
the conditional probability of finding, at a scale Q 2 , a parton i with momentum fraction x in 
a proton P while a hadron h with momentum fraction z is detected. Fracture functions obey 
non-homogeneous evolution equations [9]: 
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The first term in the above equation (see Fig.(3a)) describes the evolution of the active parton j 
while the hadron h is detected. The second term (see Fig.(3b)) takes into account the production 
of a hadron h by a time-like cascade initiated by the active parton j. The Pj' l (u ) represent the 
unsubtracted Altarelli-Parisi splitting functions [24], Since perturbative evolution is at work even 
in target fragmentation region, we expect that a non negligible amount of transverse momentum 
is produced there. We thus generalize these distributions to contain also transverse degrees of 
freedom. The fracture functions M. Ph {x, k p , z,p±,Q 2 ) give the conditional probability to find 
in a proton P, at a scale Q 2 , a parton with momentum fraction x and transverse momentum k± 
while a hadron h, with momentum fraction z and transverse momentum pj_, is detected. Under 
these assumptions the following evolution equations can thus be derived [27]: 
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As in the longitudinal case two terms contribute to the evolution of TMD fracture functions as 
displayed in Fig. (3). The homogeneous one has a pure non-perturbative nature since involves 
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the fragmentation of the proton remnants into the hadron h. The inhomogeneous one takes into 
account the production of the hadron h from a time-like cascade of parton j and thus is dubbed 
perturbative. The TMD fracture functions fulfil the normalization condition: 




Figure 3: Evolution of fracture functions M: a) homogeneous term; b) inhomogeneous one. 


J d 2 k ± J d 2 pj_Mp h (x, k±, z,p±, Q 2 ) = Mp h (x, z, Q 2 ), (13) 

as direct consequence of the kinematics of both terms in the evolution equations, eq. fll 2jl . The 
proof of factorization, i.e. that all singularities occuring in the target remnant direction can be 
properly renormalized by the less inclusive TMD quantity M.p h (x, k±, z, J>j_, Q 2 ), is still lacking 
at present. In the following we assume such a factorization to hold. 


IV. Solutions 

TMD evolution equations © and ® can only be approximately diagonalized by the joint 
Fourier-Mcllin transform 

V n (b,Q 2 ) = J d 2 p ± J e ^dz z n V(z,Q 2 ,p ± ), (14) 

where b is the transverse momentum Fourier-coniugated variable [15]. An exact diagonalization is 
prevented by the kinematics structure of the distribution under integral in eq. © since it combines 
longitudinal momentum fractions with transverse momenta. Such an exact diagonalization can 
be however obtained in the soft limit, i.e. when the variables x or z approach the edge of phase 
space. In this case for the non-singlet time-like unintegrated fragmentation functions the solution 
reads [12] : 

V(z, Q 2 , Pl _) = V(z, Q 2 ) Q{Ql Q 2 , z, p ± ), (15) 

where the scale Ql sets the upper limit of the non-perturbative regime. We recall that in this 
limit the convergence of the perturbative series can be further improved by taking into account 
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soft gluon radiation enhancements. The form factor Q can therefore be computed to leading 
logarithmic accuracy [13-15] by simply demanding that a s (Q 2 ) —* a s (p 2 ± ). The expression to 
next-to-leading logarithmic accuracy has been given in Ref. [19]. These issues have been recently 
specialized to the case of SIDIS processes in the current fragmentation region in Ref. [20]. 
Away from the soft limit, the factorized structure of the solution, ea. dlOl) . is not automatically 
preserved. In this case numerical methods have shown to be useful in order to solve the equations 
in the x, z G 0(1 ) range [17,27] and to extend the solutions to be valid in the flavour mixing 
sector. As in the longitudinal case, distributions at a scale Q 2 > Ql are known if we provide a 
non-perturbative input density at some arbitrary scale Q g. We assume as initial condition the 
usual longitudinal fragmentation distribution times a ^-independent, flavour independent factor: 


O>i(z,Q 2 0 ,P±) 


V?(z,Q 2 0 ) 


q-P 2 l/<P 2 l> 

7 T < p\> 


(16) 


The gaussian p^-distribution is used to model partons intrinsic momenta inside hadrons [21], 
This issue has also been considered in Refs. [22,23]. We now discuss the phenomenological 
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Figure 4: a) semi-inclusive DIS l + P —> l + h + X process; b) kinematics of the reaction in the 
photon-proton center of mass frame. 


implications of these results. According to Ref. [18], single hadron cross-sections are usually 
parametrized in terms of four independent structure functions H i= x j _^(xB,Zh,Q 2 ,Ph±). The 
hadronic invariants of SIDIS processes are constructed by using external particles four-momenta, 
as displayed in Fig. (4a), and therefore the structure functions Hi are supposed to describe both 
current and target hadron fragmentation mechanisms. 

Once TMD evolution equations are solved, these predictions can be compared with semi- 
inclusive DIS data, as for the longitudinal case, provided that a factorization theorem holds even 
for TMD distributions. Such a theorem has been shown to hold in the current fragmentation 
region for the structure function H 2 in Ref. [16] : 

H 2 (xB,z h ,P h ±,Q 2 ) = ^ e l [ d 2 k ± d 2 p ± 5 {2) (z h k ± + p ± - P hX ) ■ 

i=q, q 

' Pp(xBi h-Fi ^J-7 ) Pi{ z h, 0(Q 2 , /i^, p 2 D ) , (17) 
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where the standard SIDIS variables are defined as: 


Zh 


PP h 
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and fi F and p 2 D are the factorization scales. The above results are accurate up to powers in 
(P 2 ± /Q 2 ) n for soft transverse momenta Ph± — A qcd- Evolution equations for T and T> are given 
in eqs. 0 and The factor C is the process-dependent hard coefficient function computable in 
perturbative QCD and to LLA we can set C=l. Provided that factorization holds for the TMD 
fracture functions, we may add, according to eq. m, their contributions to H 2 : 


H 2 {x b , z h , P h ±, Q 2 ) = ^2 f d 2 k±d 2 p_ L ^S 2 (z h k_ L +P±~ Phi. ) ■ (19) 

i=q,q 

■P l p(x B , Q 2 , k ± ) T>t(z hl Q 2 ,P±)A( 0) + (1 - x B )M t Pjh ( L x B ,k±,z,p±,Q 2 )8 2 (p± - P h± )A( l)j 


where we have identified all the three factorization scales with the hard scale, Q 2 = p? F = 
p 2 D = ji 2 M . Although formally the two contributions are simply added in eci. (fTUl) . at LLA and 
in photon-proton center of mass frame, the produced hadrons are mainly distributed in two 
opposite emispheres. Target fragmentated hadrons are produced mainly in the 9 = n direction 
while current fragmented hadrons mainly along 9 = 0 direction, ffere 9 is the angle of the 
produced hadron h with respect to the photon direction, as shown in Fig. (4b). In order to 
keep track of the emission angle of the detected hadron h, we supplement current and target 
Lamentation terms in eq. m with an angular distribution A(v) [ 11 ]. The angular and energy 
variables v and z are defined as: 

E h 1 — cos 9 

E p {i-x B y 1 2 

In eq. 03 ), Eh and E p denote respectively the energies of the detected hadron and of the incoming 
proton in the photon-proton center of mass frame. The variables 0 and v are a useful frame- 
dependent representation for the hadronic invariant Zh in two respects: z reduces to Zh in the 
current fragmentation region so that we recover the standard definitions, while for low ^-values 
we can distinguish soft hadrons (z —> 0) from the ones produced in the target remnant direction 
(9 — > 7r). Since to LLA all sources of transverse momenta contributing to Ph _l have been taken 
into account we may pictorially represent eq. m as in Fig. (5). 



V. Conclusions and perspectives 

In this work we have extended ordinary distributions to include transverse degrees of freedom 
both in the current and in the target fragmentation region of semi-inclusive DIS. As long as a 
factorization theorem holds for transverse momentum dependent fracture functions, the semi- 
inclusive cross-sections are thus predictable on the whole phase space of the detected hadrons. 
Although this extension may have its own theoretical relevance, on the phenomenological side it 
also improves our knowledge of both the perturbative and non-perturbative dynamics of partons. 







Figure 5: Sources of transverse momentum in the current (left) and in the target (right) fragmentation 
region in SIDIS. Dark blobs symbolize hard partons emission. Transverse momentum Pf L ± of the detected 
hadron h is also indicated. F. D and M represent parton distribution, fragmentation and fracture 
functions respectively. 


Evolution equations and m can be straightforwardly extended also to the polarized case. 

In the target fragmentation region the evolution equations for longitudinal polarized fracture 
functions have been derived in Ref. [25]. We stress that the accuracy of eqs.©, © and m is 
set by the accuracy of the Altarelli-Parisi splitting functions. Splitting functions for partons and 
fragmentation functions are well known. They have been recently calculated in the target region 
at two loop level in Ref. [26]. Therefore the evolution equations in the current and in the target 
fragmentation region can be set as to allow the analysis of the semi-inclusive cross sections to 
the same level of accuracy. Detailed derivation of the results presented in this work and of the 
phenomenological implications are postponed to a forthcoming paper [27]. 
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